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Abstract 

This paper concerns the characterisation of second order marginals for random sets 
in a discrete setting. Under the instance of unit covariances, this problem possesses a 
combinatorial symmetry, exploited jointly in the companion paper to give a heuristic 
procedure to check realisability. In particular we disprove Matheron's conjecture, and 
explicit partially the structure of the convex body formed by realisable unit covariances 
in a finite set. 
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1 Introduction 

Let A 1 be a set, and call = {—1, 1} X the class of binary functions on X. Call unit field 
on X a random element X in The unit covariance of X is then denned by 

ply = BX x X y = 2P(X X = X y ) - 1, x, y E X. (1.1) 

The central question here is the inverse realisability problem, given a bivariate function 
{px,y'i x, y G X) in the class &x of symmetric functions on X, to check whether it can be 
realised by some unit field X (i.e. p = p x ). The typical example of a non-realisable, or 
non-admissible, symmetric function p is, with X = {1,2,3}, 

Pi,2 = 1, Pi,3 = 1, P2,3 = -1, Pi,i = 1;* = 1,2,3, 

because if p = p x for some unit field X, then X± = X 2 a.s., X 2 = X 3 a.s., but X± = — X 3 
a.s.. 

Measurability issues don't matter in this paper, so consider that B^- is endowed with the 
discrete topology. The affine transformation X i— > 2X — 1 transfers unit fields to {0, l}-valued 
processes, assimilable to random sets. Such marginal problems arise in many contexts, the 
article [3], for instance, calls contextual such an admissible function, in the sense that it can 
be inserted in a real physical context, with applications in the fields of information theory, 
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game theory, quantum mechanics. The authors use an entropic approach to solve some par- 
ticular related questions. Our goals turn more towards materials science and geostatistics, 
where such a characterisation could serve many purposes in modelisation and estimation, 
see the companion paper [9] for more details. On a more theoretical level, characterising 
the class of second order characteristics would provide insights on a spectral theory for ran- 
dom sets, which apart from Koch et al. [8] is currently a gap in the literature. Matheron 
identified a combinatorially compact description of the class of realisable unit covariances if 
card(X) < 6, highlighting an elegant combinatorial structure. 

Call tftx the class of admissible unit covariances on X x X. It is clear under the form (11. ip 
that to belong to ^x a function p has to be semi-definite positive, and lie in the space &\ of 
functions taking the value 1 on the diagonal (p x ,x = 1 for all x G X), but those conditions are 
not sufficient. Let &' x be the space of symmetric functions with finite support, and define 
for ae J? x 

g„(u) = a x>y u x u y ; u G 

x,yeX 

A necessary and sufficient condition for p G &x to be a unit covariance is that for every 
function a G 

J; Otx,yPx,y > : = Ulf g«(u), (1.2) 

z — ' ueB^. 
x,yeX 

because then the operator that p induces on the vector space spanned by constant functions 
and functions {g a : a G &'x) is positive, and one can apply Kantorovitch and Riesz-Markov 
theorems (see [ID] ) (the necessity of (II. 2p is a straightforward consequence of the positivity 
of the mathematical expectation). In this context, p is sometimes loosely referred to as being 
positive. 

The inequalities (II. 2p are the linear inequalities determining the convex set %xi and we 
show in this article how this approach, effective in [TSj and [T7], can be fruitful. Even if 
X is finite, the positivity conditions (11.2j) involve a priori an infinity of relations, and it is 
not even known what is the right-hand side of (II. 2p for general a. In the discrete setting 
and in a slightly different formulation, McMillan [T5J has proved that it is enough to satisfy 
(ll.2p for a corner positive, and x a = 0, but the corner positivity is poorly understood; it 
has been numerically studied in [T7] if X has less than 7 elements. At the present time it 
is not possible to give an algorithmic procedure that determines the positivity of a given 
function (p x ,y), and it seems like a very difficult challenge with the tools available, therefore 
the consensus is to give necessary conditions as sharp as possible for the admissibility of p. 
This problem has been posed by MacMillan [15] in the field of telecommunications. It is 
more or less implicit in many articles, and has been to the author's knowledge first addressed 
directly by Shepp [TS], and more recently by Quintanilla [17J. A series of works by Torquato 
and his coauthors (see [7j, and [19] Sec. 2.2 and references therein), in the field of materials 
science, gather known necessary conditions and illustrate them in many 2D and 3D models. 
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This question was developed alongside in the field of geostatistics; Matheron [TJ] has found 
via arithmetic considerations a wide class of necessary conditions, that he has proven to be 
sufficient if X has cardinality less or equal to 5, and he has conjectured these conditions to 
be sufficient for any X . Those conditions form the widest known class of necessary condi- 
tions used in related practical problems. Some other authors do not attack frontally this 
question, but address the realisability problem within some particular classes of models, e.g. 
Gaussian, mosaic, or boolean model (see [TJ [21 [HI EES])- In the companion paper [9], we use 
the theoretical results of the present work to give a heuristic algorithm allowing to discard 
some inadmissible covariances, and other applications. 

Our method is based on a direct study of which is a convex polytope if X has 
finite cardinality, and can therefore be applied tools from convex geometry. In Section 11.11 
we provide some structural information on ^x-, such as its dimension and extreme points, 
and from there on X is essentially assumed to have finite cardinality. As an application 
we give a bound for the number of different states of a random field realising an arbitrary 
admissible covariance. Section [2] revolves around the linear inequalities characterising ^x^ 
i.e. finding the hyperplanes supporting its facets if X has finite cardinality; we give new 
necessary conditions that are the support of the heuristic algorithm developed in [9j, and 
disprove numerically the conjecture of Matheron on the form of the supporting hyperplanes 
normals. Section [3] is devoted to more theoretical facts about ^x-, such as the elements of 
its boundary and its graph structure for faces of low dimensions. 

1.1 Realisability and convex geometry 

It is a fairly trivial fact that ^i/x is convex, indeed if X\ and X2 are two unit fields of X, 
then the segment [p Xl , p X2 } is comprised in %x because for every t G [0, 1], tp Xl + (1 — i)p X2 
is the unit covariance of 



where B is an independent Bernoulli variable with parameter t. The convexity is exploited 
here to lay out some notation and basic facts about realisability. See the appendix at the 
end of the paper for notation and basic vocabulary in convex geometry. For a function v on 
X, denote (v £g> v) x ,y = Vx^y, v <E> v G If X is a unit field, then p x = E[X <g> X], where 
the expectation is taken component-wise. Furthermore the structure of vertices of ^x can 
be made explicit. 

Proposition 1.1. The extreme points of °}/x are the u <8> u, u G Bx- 

Proof. Take u G B^, and assume that u <E> u = tp Xl + (1 — t)p X2 for some unit fields Xi, X 2 
and t G [0,1], t ^ 0,1. Since for every x,y G X we have tp Xl y + (1 - t)p X2 y G {-1,1}, 
necessarily p x y = p x y G {—1, 1} for every pair (x,y) of X 2 . It follows that both Xi and X 2 
are deterministic with p Xl = p X2 , and therefore u <g> u is an extreme point. 
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Conversely, take p = p x G 'Wx, assume that for some x, y G X, t = p x y G (—1,1), 
meaning P(X X ^ X y ) ^ {0, 1}, and define the two unit fields Xi and X 2 with respective 
laws 

P(Xi G A) = P(X G A | X, ^ Xj,) 
P(X 2 G A) = P(X G A | X x = X y ) 

for A a subset of B^>. We have, by conditioning, 

p x =tp Xl + (1 -t)p x \ 

It follows that p x lies in the relative interior of the segment [p Xl , p X2 } (this segment is not a 
singleton because 1 = p Xl y ^ = 0). Thus p x is not an extreme point of tyfx- D 

In the sequel we focus on the finite case X = [N] = {1,2,..., N}, N > 1. For the sake of 
clarity index [N] is replaced with N. Functions of J^/v are identified with N x N matrices. 

Remark 1.2 (Notation for triangular arrays). Call 7r the projection operator that takes 
a matrix p G J^jv to its supra-diagonal components 7r(p) = (p^; 1 < i < j < N), and 
call ^ N = 7r(J^v) the space of supra- diagonal triangular arrays. An element p of & N is 
represented as the triangular array 

/ Pl,2 ■ ■ ■ Pi,n ^ 

Pn-i.n 

\ J 

Since any p G automatically has the diagonal filled with l's, it is somehow intuitive 
that mapping %y to = 7r(^jv) does not lose any relevant information. The dimension 
of ^ N is (In '■= N 2 > anc ^ we P rove below that % N is a full- dimensional convex subset of 
The elements said to be under the product form are those that can be written, for some 
v G R N , 

v ® v := 7r(v <g> v) = (viVj)i<i<j<N- 

Proposition 1.3. ffl N has a non-empty interior in & N . Its extreme points are the 2 N ~ l 
vertices of the form u Q u for u G B^r. 

Proof. It suffices to prove that ^ N contains at least linearly independent vectors. Take 
i < jo in [N]. Let X k , k ^ j be independent Rademacher variables (P(X = 1) = PpQ = 
— 1) = 1/2) and define the unit field X = (X 1; . . . , Xj Q _i, X io , Xj 0+1 , . . . , X N ). Its unit 
covariance is the canonical vector e io j = (l{ i=i0)j=jo })i<i<j<N- Therefore & N contains and 
the djsr such canonical vectors, and has nonempty interior in 

In virtue of Proposition [TTTJthe extreme points of %v are the u® u for u G Bjy, whence the 
extreme points of % N = 7r(%v) are of the form u Q u, u G Bjy. Conversely we can prove that 
each uQu is an extreme point by mimicking the arguments from the proof of Proposition II. 1[ 
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simply by assuming x < y. All the uSu are hence extreme points, and we have u u = v v 
if and only if u = v or u = —v. It follows that the number of extreme points is the number 
of binary vectors, under the identification u = — u; we arrive at 2 7V ~ 1 . □ 

A consequence of these remarks is that the number of different values taken by a random 
set realising a given unit covariance can be chosen to be no larger than + 1. 

Proposition 1.4. Any unit covariance p £ %v can be realised by a random field X that 
takes at most + 1 distinct values. 

Proof. Applying the Minkowski-Caratheodory theorem to in Tx(p) can be expressed 
as the convex combination of d^ + 1 extreme points u k u k , u k £ B N , 1 < k < d N + 1 

d N +l 

n (p) = Yl pkuk ufe ' 
*!=i 

with pk > 0,YlkPk = 1; meaning p is the unit covariance of the unit field X which law is 
defined by 

P(X = u fc ) = Pk £ [0, 1] 
for every 1 < k < djv + 1. □ 

Remark 1.5. Another consequence is that contains a non-empty open ball. Using 
Gaussian unit covariances, it is proved in [9] that the euclidean ball of centred in with 
radius \/2/ir is comprised in 'W^. 



2 Checking realisability numerically and Matheron's 
conjecture 

This section focuses on the practical problem of checking, when X = [N],N > 1, is finite, 
the validity of a given function p. We saw in Proposition 11.31 that the problem is properly 
posed in J^, and since is a polytope, it can be written as a finite intersection of half 
spaces, 

h N 
k=l 

where is the number of facets, the p k £ are the outer normals of the Ok £ K are 
the corresponding offset values, and 

W;t = {pe^:(p,» k )<o k }. 

For a finite-dimensional vector space V, call (■, - )y the canonical scalar product on V. All 
it takes to find Ok knowing p k is to compute the infimum of (p k , -)jr* on the 2 N ~ l extreme 
points of (prop. II. 3p . The computational problem consists in finding the normals p. , 
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the complexity of computing the corresponding offsets is much smaller, see Tabled! In [14j, 
Matheron studied the structure of the convex set constituted by realisable covariograms, 
similar to that of For iV < 5, he computed by hand the normal vectors and detected a 
recurring pattern; he then conjectured that this form should be valid for every N > 1. We 
transposed and extended his results with numeric computations here to unit covariances. 
Introduce the class 

N 

Etv = {e G Z N : Ujej = 1 for some binary vector u G B N = { — 1, 1}^}. 

i=l 



The following theorem gives a neat characterisation of unit covariances for N < 6. 

Theorem 2.1. (i) Take N < 6. Every outer normal p k ,l < k < h^, can be written 
under the form e k <S> e fc for some e G E^. A matrix p G belongs to %v if and only 
if for every e G Ejy, 

N 

(ii) Let N > 1, and p G For every e G Z N with odd sum, p satisfies 112.1}) . 

Proof, (i) The outer normals p k of can be listed by a linear programming algorithm 
which input is the 2 N ~ l extreme points of and output is the outer normals and corre- 
sponding offset values of the polytope spanned by these extreme points. We used the cdd+ 
algorithm [I], the list of these normals can be found at 

http://www.math-info.univ-paris5.fr/ rlachiez/realisability. We then checked 
that each normal p k could indeed be put under the form e fc e k for some e k G E^. Thus 
p G J^jy is realisable if and only if for every 1 < k < 

(7r(p), e k S) e fc ),^ > inf (•, e k e fc )^. (2.2) 

Since the minimum of a linear form on a convex polytope is necessarily reached in an extreme 
point, Prop. 11.31 yields that (12. 2 p is equivalent to 

(tt(p), e k e k ),?* > inf (u u, e k e fc )jr* , 1 < k < h N . 

N ugB]v N 

Using the fact that p G J^, £) i=1 Pu(e k ) 2 = Y,i=i UjU £ (ef) 2 for every u G B^, 1 < k < h N , 
the prior condition is therefore equivalent to 

(p, e k e k }^ N > inf (u u, e k e fc )^ = inf ( V Ui eA = 1 
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because J2i=i u i e i * s an °dd number, and the value 1 is reached for the particular j e B w 
such that Y2i=i e i u i = 1- Finally, a matrix p G indeed belongs to if and only if 

JV 

Pije^ej > 1 for every 1 < k < h N . 
The necessity is contained in point (ii). 

(ii) If p = px for some unit field X, and e G Z N is such that £\ is an odd number, 

(px, e ® e),^ = E ^ > 1 

because £\ e^Xj is a.s. an odd integer. □ 

Remark 2.2. 1. Even though Ejv is infinite, point (i) provides a finite time procedure to 
check the realisability of a given matrix p, see the details in the companion paper [9]. 

2. Any vector of E N has odd sum, thus point (ii) indeed applies to E N . 

3. Matheron only made computations up to N = 5, for which he found that we can 
furthermore state that G { — 1,0,1}, but numerical computations showed that this is 
no longer true for N = 6. 

4. The point (ii) also originates in |14j . 

5. The point (i) concerns the normal vectors to all the facets of the polytope fy^- Math- 
eron originally formulated his result only for the normal vectors of the facets touching 
the vertex (1, . . . , 1). In the meantime he mentioned the idea that this could be used 
for other vertices by exploiting the combinatorial symmetry of unit covariances. 

We give in table [1] the number of facets for N < 7. It explodes as the dimension 



N 


dimension 


h N 


3 


3 


4 


4 


6 


16 


5 


10 


56 


6 


15 


368 


7 


21 


116 764 



Table 1: number of facets of %Z 



increases. For N = 8, the algorithm was stopped in dimension d$ = 28 after that 2 216 100 
outer normals were found, no more memory space being available. To go further in the 
computations, one has to go deeper in understanding the structure of ty^. Section [3] gathers 
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some metric and topological facts about that can help optimise the processing time and 
understand the structure of 

Theorem 12.11 provides necessary conditions that are used in the companion paper jH] to 
design an algorithm able to discard some spherical variograms as admissible covariograms, 
a recurrent problem in geostatistics (see [IT] . Sec. 3.2.2). 

2.1 Matheron's conjecture is not true 

Matheron's conjecture is appealing, as if it were true it would provide a procedure to effi- 
ciently check the realisability of a matrix for any N >1. We show in this section via theoretic 
and numeric arguments that the conjecture fails at N = 7. Recall that convexity-related 
notation and vocabulary are introduced in the appendix. 

The conjecture was originally stated for the covariogram of a random set Y C [N] 



where ly is the indicator function of Y. Calling Yn the class of admissible covariograms 
(convex for the same reasons than %v), and Y^ = tt(Yn), Matheron conjectured that the 
outer normals of cone(0; Yff) are of the form eSe for e in 7L N with Yli=i e « = 1- Now 
$(7y) := 1 — A'jy is the unit covariance of the unit field X = 2Y — 1, therefore easy 
computations show that the outer normals of cone(l; (where 1 = (1, . . . , 1) = $(0)) are 
those same eSe. We prove below that the conjecture, even under a slightlier general form, 
is equivalent to Th. 12. 1L 

Remark 2.3. Matheron conjectured that conditions (12.11) characterise matrices contained 
in cone(0; Yfi), more precisely that for 7 G J^ N , 



for every e e Ejv, but not that they characterise the whole convex Y£ as can be seen in the 
literature. 

Theorem 2.4. For any N > 1, Matheron's conjecture is equivalent to the conjecture that a 
matrix p G 1S a un ^ covariance if and only if it satishes ( 12. Jj) . 

Proof. Let 7 G = {7 G ^ N : 7^ = for all i G [N]} and p G be linked up by 
p = 1 — 47, hence p is realisable as a unit covariance iff 7 is realisable as an indicator 
covariogram. 

Let us assume that Matheron's conjecture is true, and that p satisfies (12.11) ; we must 
prove that p is realisable. In particular, for e G with unit sum (i.e. J2i e « = 1) 




\P(1 Y (X) ? l Y (y)) 



i,je[N) 



7 G cone(0; Yn) if and only if Jij e i e j ^ 

l<i<j<N 



(2.3) 




l<i,j<N 



l<i,j<N 



s 



whence by Matheron's conjecture tt(j) lies in cone(0; ^). Applying $ yields that n(p) is in 
cone(l; Let us call, for u G {—1, 1}^, 9 U the transformation of J^jv defined by 

9 U fl = (UjUj//y)l<jJ<JV, p G ^jv- 

The rotated matrix 9 u p also verifies 

^e^V)^- > 1 

ij 

for e G Z w with unit sum because p satisfies (12.11) . It follows that n(8 u p) also belongs to 
cone(l; ^^), and therefore ir(p) G cone(w St/; < ^). Thus in each extreme point u Q u of tyfi, 
u G { — 1,1}^ (see Prop. II. ip . n(p) is contained in cone(u <S> u;^^). Thus n(p) is in 
meaning it is realisable. 

In the other direction, if (12.11) characterises the realisability of unit covariances, then 
cone(l; is characterised by inequalities in (12. ip that become equalities if pij = 1, i.e. 
inequalities such that £V e$ = 1. It follows by applying $ _1 that cone(0; "V£) * s characterised 
by (12. 3p . which is exactly Matheron's conjecture. 

□ 

In the light of the previous proposition, Matheron hence conjectured that every normal 
vector of is under a product form v0v. We give in the file rho7 . ine the normals of 
computed with cdd+. The following array r G is one of these normals. 

/ —2 1 1 1 1 \ 
-2 1 -2 -2 -3 
-1112 
t= -1-1 -2 

1 2 
2 

V / 

It is clear that r is not under the product form v Q v for some v G M N because T14 = 0, 
but neither its first line nor its 4-th column are filled with 0's. In conclusion Matheron's 
conjecture is not true for N = 7. We proved in a similar way that it was not true for N = 8. 

3 The structure of 

The task of giving a tractable characterisation of realisable covariances, or unit covariances, 
seems a very hard one, not to say impossible. This problem, combinatorial in nature, relies 
heavily on arguments from convex geometry to analyse We give in this section topo- 
logical facts about in order to better apprehend its structure both for a geometric and 
a graph-theoretic description. The convex body bears some peculiar properties, which 
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might lead to think that the combinatorial problem of characterising is better appre- 
hended with this geometric approach, and stands on its own as an interesting theoretical 
problem. 

With the current knowledge of one should more rely on numeric considerations to 
test the validity of a given model. Unfortunately, as is apparent in Table [H the complexity of 
such a task explodes as the cardinality iV of X increases. Understanding better the structure 
of can enable programmers to design a more adapted linear programming algorithm to 
find the outer normals of which amounts to find the necessary and sufficient conditions 
of realisability. 

Let us start by recalling that is a polytope of which all 2 N ~ 1 extreme points are 
vertices of the hypercube (Prop. 11.31 ). For algorithms based on ray-shooting queries it is 
interesting to have an interior point of or of its dual. Here plays perfectly this role, 
and this can be quantified. 

Proposition 3.1. Calling Bo(r) the euclidean ball of with radius r centred in 0, we 
have 

B (V2/n) CW*C B (^). 

Proof. The first inclusion translates Remark 11.51 itself relying on classical results about 
Gaussian covariances and exploited in [9]. The second traduces the fact that the extreme 
points of are vertices of the hypercube. □ 

The next statement locates the singular admissible covariances on the boundary of fy^. 

Theorem 3.2. (i)The unit co-variance p of a unit field X is singular if and only if for some 
non-trivial real deterministic coefficients Ai, . . . , Ajv the components of X satisfy 

^ AfeX fc = Oa.s.. 

k 

(ii) Singular unit covariances lie on the boundary of 

Proof, (i) Let X be a unit field and p its unit covariance. As a semi-definite positive matrix, 
the singularity of p is equivalent to the existence of a non-trivial family of scalar numbers 
{X k : 1 < k < N} such that 

AfcAjPfcj = 0, 

l<k,j<N 

which exactly means that E(^^L 1 X^X^) 2 = 0. 

(ii) Let p be a singular unit covariance, and A&, k = 1, . . . , N, non-trivial coefficients such 
that 

X k\p kj = 0, l<j<N. 

l<k,j<N 
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Two at least of the A& are non zero (because p has l's on the diagonal), say A^ and Aj 2 , and 
put a = sign(A il Aj 2 ) G { — 1, 1}. Let us take e > and 

Pij = Pj - a£ if {hj} = {h, ^2}, 
Pij = Pij otherwise. 
Then if p e was realisable by some unit field X s , we would have 

N N N 

E(^X k Xl) 2 = X k\p e u = Pki^kk - 2AiA 2 (X£ = -X n X V2 ae < 0, 

fc=l k,l=l k,l=l 

which is impossible. Thus there are unrealisable covariances arbitrarily close from p, it 
follows that p is on the boundary of ^ 

Remark 3.3. It is conversely not true that every p on the boundary is singular (otherwise 
finding the outer normals of and thus characterising realisable unit covariances, would 
be easy). According to Th. 13.21 it would mean that £\ \Xi = a.s. for some non-trivial 
family (A*), but if for instance X takes the values (—1, 1, 1), (—1, —1, 1), (1, —1, 1) each with 
probability 1/3 yields the non-singular unit covariance 

/ 1 -1/3 -1/3 \ 
Px = -1/3 1 -1/3 . 
\ -1/3 -1/3 1 / 

On the other hand 7r{px) lies on the border of because it is a convex combination of 
three out of the 4 vertices of the tetrahedron (see Prop. 11.31 or [H] p. 110). 

A manner of describing exhaustively the topological structure of is to explicit the 
hypergraph structure of its vertices. The following theorem states that for k < N, the k-th 
order hypergraph structure of is complete, in the sense that the simplex formed by any 
fc-tuple of extreme points is contained in ^^'s boundary. 

Theorem 3.4. For every k < N, and k-tuple u 1 , . . . , u k of { — 1, 1}^, the k-dimensional sim- 
plex formed by their respective unit covariances p u i , . . . , p u k only contains singular matrices 
and thus lies on the boundary of fyf^- 

Proof. A covariance of the form p = Yl^iPiP^ is the covariance of the unit field X with law 
P(X = u*) = Pi, allowed to take values among the k vectors u\ Since k < N, there exists a 
non-trivial family Ai, . . . , Ajv such that 

N 

for all i, whence a.s. Ylj = 0. It follows by Th. I3.2( i) that p is singular. The conclusion 
comes by applying Th. 13.21 □ 

In some sense has the topological structure of the c?Ar-dimensional simplex as long as 
one only looks at dimensions strictly smaller than N . Higher dimensional facets are harder 
to explicit. 
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Appendix: Convex geometry 

Let C be a convex set of a vector space V . The extreme points of C are the x G V such 
that, for every y, z G C, if x lies in the segment [y, z] spanned by y and z, then either x = y 
or x = z. Denote them by ext(C). All the subsequent statements are only made for finite 
dimensional space V = W 1 . If X is a subset of V, the convex hull of X is the smallest closed 
convex set of V containing X. The relative affine space span(C) of a convex C is the vector 
space spanned by C. In the topology of span(C), C has a non-empty interior and the relative 
dimension of C is defined as the dimension of span(C). Say that C is fully dimensional if 
span(C) = lR d , equivalent to the fact that C has a non-empty interior in IR d . Given a point 
x G C , denote by cone(a;; C) the smallest affine convex convex cone with vertex x containing 
C. 

Polytopes 

A polytope P of M. d is the convex hull of a finite number of points x\,. . .,x q of R d , denoted 
P = conv(xi, . . . , x q ). For < k < d, a k- dimensional face of P is any set S C dP such that 

5 is the convex intersection of P with a affine hyperplane and the relative dimension of S is 
exactly k. For any < k < d— 1, the data of all its /c-dimensional facets characterise P. Its 
0-facets are its extreme points, also called vertices, and the corresponding description of the 
form ext(P) = {x±, . . . ,x q } is called V-description. At the opposite, its (d — l)-dimensional 
facets are sometimes just called facets, and yield the so-called H- description. The latter can 
be understood via a finite family of affine half-spaces H\, . . . ,H q where 



for some G M. d called outer normal and corresponding offset value Oi, 1 < % < q. Given 
any x G ext(P), the convex cone of C generated by x is obtained by retaining only half 
spaces with x in their boundary 



For k < d, call fc-th order hypergraph of C the class of fc-tuples of extreme points of C 
{x\, . . . , Xk} such that the simplex formed by x±, . . . , Xk lies on the boundary of C. 



Hi = {yeR d : (y,tu)<Oi} 




i:{x,M}=Oi 
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